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Abstract. We study exclusion processes on the integer lattice in which particles 
change their velocities due to stickiness. Specifically, whenever two or more particles 
occupy adjacent sites, they stick together for an extended period of time and the 
entire particle system is slowed down until the "collision" is resolved. We show 
that under diffusive scaling of space and time such processes converge to what one 
might refer to as a sticky reflected Brownian motion in the wedge. The latter 
behaves as a Brownian motion with constant drift vector and diffusion matrix in 
the interior of the wedge, and reflects at the boundary of the wedge after spending 
an instant of time there. In particular, this leads to a natural multidimensional 
generalization of sticky Brownian motion on the half-line, which is of interest in 
both qucueing theory and stochastic portfolio theory. For instance, this can model 
a market, which experiences a slowdown due to a major event (such as a court trial 
between some of the largest firms in the market) deciding about the new market 
leader. 



1. Introduction 

Stochastic processes with sticky points in the Markov process sense have been 
studied over the course of the last three decades. In the now classical papers [8] 
and pp, the authors analyze sticky Brownian motion on the half-line, which is the 
process evolving as a standard Brownian motion away from zero and reflecting at zero 
after spending an instant of time there (as opposed to a reflecting Brownian motion, 
which reflects instantaneously at zero). These papers show that sticky Brownian 
motion arises as a time change of a reflecting Brownian motion, and that it describes 
the scaling limit of random walks on the natural numbers whose jump rate at zero 
is significantly smaller than the jump rates at positive sites. 

In stochastic analysis the stochastic differential equation (SDE) 

(1-1) dS(t) = l { s W >o} dB{t) + rj l {S (t)=o} dt 

satisfied by sticky Brownian motion has drawn much attention, as it is an example 
of a SDE for which weak existence and uniqueness hold, but strong existence and 
pathwise uniqueness fail (see [5]). In fact, in [18] (see also the survey [6]) it is shown 
that a weak solution to (11. ip cannot be adapted to a cozy filtration, that is, a filtration 
generated by a finite or infinite-dimensional Brownian motion. 
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The present study is motivated by the question of how one can define and analyze 
multidimensional analogues of (11 .ip and whether solutions to the corresponding sys- 
tems of SDEs arise as suitable scaling limits of interacting particle systems in analogy 
to the findings of [8] in the one-dimensional case. In [TTj Section 3] it is shown that 
a large class of reflecting Brownian motions in the n- dimensional wedge 

W = {x eW 1 : xi < x 2 < . . . < x n } 

arise as limits of certain exclusion processes with speed change under diffusive rescal- 
ing. As we show below, sticky Brownian motions with state space W can also be 
obtained as scaling limits of suitable exclusion processes with speed change. 

1.1. Exclusion processes with sticky particles. To simplify the exposition, we 
next describe a simple class of particle systems which converge to sticky Brownian 
motions in W in the scaling limit, and postpone the description of the much wider 
class of particle systems that we can handle to section [31 We fix the number of par- 
ticles n G N, and also rate parameters a > 0, <d L = {&i,j) ie ^ j- e r n _i] e IP> oo) nX ( n_1 \ 
and 6^ = (0?,). . . . . G [0, oo)" x(n ^ 1) with the notation 

[n] = {l,2,...,n}. 

For a fixed value of the scaling parameter M > 0, the particles move on the rescaled 
lattice Z/ vM ; to describe their motion we introduce the following Poisson processes, 
all of which are independent, and all of which have jump size h=. For % G [n], the 
Poisson processes Pi and Qi have jump rates Ma, while for i G [n], j G [n — 1], the 
Poisson processes Ly and Rij have jump rates and yM9^j, respectively. 

In addition, for notational convenience we introduce ghost particles at ±oo, namely: 
X A/ (•) = -oo and X^ +1 (•) = oo. For any initial condition Xf (0) < X 2 M (0) < 
. . . < X^ 1 (0) on Z/\/M we can then define a particle system evolving on Z/a/M in 
continuous time by setting 



AX? (t) = d (P t (t) - Qi (t)) 

n- 1 

+ E 1 {xf i t )+ ^ ff<X M l{ t ) ,xf(t) + ^ M =xf +1 (t)} d ^,i (*) 

7 = 1 



n-1 



1^ 1 {x?L 1 {t)+^<x?'{t)^{t)+-^=xfl 1 {t)} dL M (*) ' 
i=i 

for i G [n]. Note that (II. 2p guarantees that for any t > 0, the particle configuration 
(X^(t),X^(t), X^(t)) is an element of the discrete wedge 

W M = Ix G (z/y/My : x k + -1= < x k+1 , k G [n - 1] J . 

Intuitively, when apart, the particles move independently on the rescaled lattice 
Z/ \/M according to the processes Pj — Qj, z = 1, 2, . . . , n (in particular, with jump 
rates of order M); however, when two particles land on adjacent sites — an event we 
describe as a "collision" — the system experiences a slowdown: the particles change 
their jump rates to the ones of the processes L it j and Rij, i G [n], j G [n — 1], which 
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are of order \/M. The interaction between adjacent particles can be described as 
stickiness, as it takes a long time (on the time scale Mi) until the collision is resolved 
and the particles return to jump rates of order M. 

1.2. Convergence to multidimensional sticky Brownian motions. The de- 
scribed particle systems converge to a sticky Brownian motion in W under the fol- 
lowing assumption. 

Assumption 1. Define the speed change matrix V = (vi t j)i £ [ n ] t j £ [ n -i] by setting 
Vij = 6fj — Qfj, and the reflection matrix Q = (qj,j')j,j'e[n-i] by setting q^y = 
Vj+ij' — Vjji. When there is a collision between particles f and f + 1 and no other 
collisions, then the velocity of particle i is given by Vij>, and the velocity of gap j 
between particles j and j + 1 is given by q^y. Assume that the matrix Q is completely- 
S, in the sense that there is a A G [0, oo)™" 1 such that Q A G (0, oo)™" 1 and the same 
property is shared by every principal submatrix of Q (see [17] for several equivalent 
definitions). 

Under Assumption [Tj, we have the following convergence result. 

Theorem 1. Suppose that Assumption [1] holds, and also that the initial conditions 
{ [X^ 4 (0) , XJ} 4 (0) , . . . , X^f (0)) , M > 0} are deterministic and converge to a limit 
(xi,X2, . . • , x n ) G W as M oo. Then the laws of the paths of the particle systems 
{ (Xf 1 (•) , X 2 M (•),•••, X^ 1 (•)) , M > 0} on D ([0, oo), R n ) (the space of cadlag paths 
with values in W 1 endowed with the topology of uniform convergence on compact sets) 
converge to the law of the unique weak solution of the system of SDEs 

n-l 

(1.3) dXi (t) = l {Xl (t)<x 2 (t)<...<x n (t)} V2adWi (t) + ^ l{Xj(t)=x j+1 (t)} dt, 

3=1 

i G [n], in W starting from (xi,x%, . . . , x n ). Here (Wi, W2, • • • , W n ) is a standard 
Brownian motion in W 1 . 

The solution to (11. 3p evolves as a Brownian motion when away from the boundary 
dW ofW, does not spend a non-empty time interval on dW, however satisfies 

P (£ ({t > : X(t) G dW}) > 0) > 0, 

where £ is the Lebesgue measure on [0, 00). 

We refer to the solution of (ll.3p with a = 1/2 as sticky Brownian motion in 
W with reflection matrix V. We choose this terminology because the SDE (11. 3p 
generalizes one-dimensional sticky Brownian motion as in [El H], and also because 
it is consistent with the terminology used in [T7] and the references therein dealing 
with instantaneously reflecting Brownian motions. 

In section [3] we prove a much stronger result than Theorem [U allowing for non- 
exponential interarrival times between the jumps in the processes Pj, Qi, L it j and Rij 
as well as for dependence between the latter processes (see Theorem [6]). This then 
leads to the definition of a sticky Brownian motion in W whose components have 
unequal drift and diffusion coefficients. In addition, it is not hard to see from the 
proof that for each jump parameter or 0fj which is zero, we can choose the jump 
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rate of the corresponding process Ljj or R it j to be of order o{\fM^j (not necessarily 
identically zero) for the result of Theorem Q] to still hold. 

One of the main technical difficulties in the proof of Theorem [T] and its extension 
(Theorem O in section [3]) is posed by the indicator function appearing in the dif- 
fusion matrix of the limiting process. This is in contrast to the main convergence 
result in [TTJ, where the martingale part of the limiting process is a Brownian motion. 
Another major difference to the setting in [11] is that we consider the full class of 
completely-5 reflection matrices and are dealing with weak solutions of the limit- 
ing stochastic differential equation; whereas in [11] only a special class of reflection 
matrices is considered, allowing for a pathwise construction of the limiting object. 
Finally, we allow for dependence of interarrival times between jumps for different 
particles in Theorem below, which has not been addressed in [TT] . 

1.3. Applications. We mention two potential areas of applications for the process 
in ( II. 3p and its extensions that appear in section [2J To this end, we recall that 
reflected Brownian motions in W give a class of tractable descriptive models for 
the logarithmic market capitalizations (that is, the logarithms of the total values 
of stocks) of firms in a large equity market (see, e.g., [H]). These models lead to 
realistic capital distribution curves in the long-run and are also able to produce a 
realistic pattern of collisions. In the same spirit, one can think of ( II. 3p as a model 
for the logarithmic market capitalizations in an equity market in which the market 
experiences a slowdown whenever there is a possibility that two firms will exchange 
their ranks (described by a collision). For example, one can imagine a court trial 
between two firms, the result of which decides which firm becomes the market leader, 
leading to a slowdown of the market right before the time of the verdict as the market 
participants await the result of the trial. 

Another area of potential applications is the study of diffusion approximations of 
storage and queueing networks. It is well-known (see, e.g., the survey [19] and the 
references therein) that reflected Brownian motions in the orthant describe the heavy 
traffic limits of many queueing networks such as open queueing networks, single class 
networks and feedforward multiclass networks. Moreover, in [5] the authors explain 
how sticky Brownian motion on the half-line can be obtained as the diffusion limit of 
modified storage processes, and in [T3] a related single server queueing system is stud- 
ied. In view of these results, we expect sticky Brownian motion in the orthant Ifjr 1 , 
given by the spacings process 

(X 2 (•) - X 1 (•) , X 3 (■) - x 2 (■),..., X n (•) - X„_! (•)) , 

to arise in the diffusion limit of suitable queueing networks as well. 

1.4. Future directions. A natural direction for future work is to study other types 
of sticky interaction between particles. Even in the class of exclusion processes in one 
dimension, there are avenues to be explored. For instance, the exclusion processes 
described by (11.21) experience a global slowdown when a collision occurs, whereas 
for some applications it would be interesting to consider particle systems with local 
slowdown. We believe that the techniques we develop in section [3] would carry over 
to such a setting with appropriate modifications; however, the difficulty of proving 
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convergence of such processes to the appropriate continuous object comes from prov- 
ing uniqueness for the limiting SDE. We expect the solution of this SDE to spend 
a positive amount of time on lower-dimensional faces of the wedge W, making the 
analysis of the process more difficult. 

1.5. Outline. The rest of the paper is structured as follows. Section [2] is devoted to 
the study of sticky Brownian motions in W. In subsection 12.11 we give the proof of 
existence and uniqueness of the weak solution to a system of SDEs generalizing (11. 3p . 
Then, in subsection 12. 21 we show that the solution is a Markov process and study the 
invariant distributions of a suitably normalized version thereof. Subsequently, sec- 
tion [3] deals with the convergence of exclusion processes to sticky Brownian motions 
in W. In section 13.11 we first prove Theorem [H and then in section 13.21 we state and 
prove our main result, namely a generalized version of Theorem [U which deals with 
the convergence of exclusion processes with non-exponential and possibly dependent 
jump interarrival times to sticky Brownian motions in W. 

2. Multidimensional sticky Brownian motions 
This section is devoted to the study of the system of SDEs 

n-l 

(2.1) dXi (t) = l {Xl (t)<x 2 (t)<...<x n (t)} {k dt + dWiit)) + ^ l{ Xj (t)=x j+1 (t)} u*j dt, 

i=i 

i G [n], where b i: i G [n] are real constants, W = (Wi, W 2 , ■ ■ ■ , W n ) is an Tri- 
dimensional Brownian motion with zero drift vector and a strictly positive definite 
diffusion matrix <£ = (c^/)^/^], and V = (vi,j)ie[ n ], je[n-i] is a matrix with real en- 
tries. We note that the diffusion matrix of the process X is both discontinuous and 
degenerate, so neither existence nor uniqueness of a weak solution to (12. ip can be 
obtained directly from the classical results in [16] or [2]. 

2.1. Existence and uniqueness. In this subsection, we show that Assumption [T] is 
necessary and sufficient for the existence and uniqueness of a weak solution to (12. ip . 
Furthermore, even under Assumption [1] one cannot expect a strong solution to exist. 

Theorem 2. Under Assumption^there exists a unique weak solution to (12. ip . More- 
over, if Assumption^ does not hold, there is no weak solution to (12. ip . 

Proof. Step 1. We start with the proof of weak existence. To this end, we first apply 
the main result in |T7] to deduce that there exists a weak solution on a suitable 
filtered probability space (J-f) t>0 , P} to the following system of SDEs: 

n-l 

dZ t (t) = (b i+1 - h) dt + dBi (t) + QiJ dA j (*) , ie[n-l], 

3=1 

where B = (Bi, B 2 , . . . , B n _i) is a Brownian motion with zero drift vector and diffu- 
sion matrix A = (at ) t/) iji / e [ n _i] given by 

(2.2) a^j' = dji + Cj+i^'+i — Cj^'+i — Cj+i^/, 

and the Aj (•), j G [n — 1] are the semimartingale local times at zero of the processes 
Zj(-), j G [n — 1], respectively. Note that here we have used the fact that the 
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matrix Q is completely-iS (see Assumption [1]). Next, we can find (after extending the 
underlying probability space if necessary) a Brownian motion = (/3i,/3 2 , . . . , 
with zero drift vector and diffusion matrix £. such that 

£»(•) = A+i(-)-A(0> ie[n-l]. 
Therefore, we can define X = (Ai, X 2 , . . . , X„) as the unique process satisfying 

n n n— 1 

]T x,(t) = ^ (m + A(t) + ]T ^,A(*)) , 

(A 2 (t) - A\ (0 , . . . , X n (t) - X n _, (t) ) = (Z 1 {t),..., Z n _ x (t) ) , 
for all t > 0. Finally, we let 

n-l 

T{t) := t + A(t) := t + ^Ajit), t>0, 

r(t) : = inf{s > : T(s) =t}, t> 0, 
and set X (•) = X (r (■)). Then clearly 

„ T (.) n-l 

X i (-)-X i (0)= / 6,df + A(r(-))+y;^A i (T(-)), ie[n}. 
Jo J=1 

Moreover, we note that r(-), A(r(-)) are nondecreasing functions, which induce 
nonnegative measures dr (•), dA (r (•)) on [0, oo) satisfying 

(2.3) dr(t) + dA(r(t)) = dt. 

Therefore, for i G [n] we have 

r(-) rr(-) 

bidt= / bi l{x 1 ( t )<x 2 (t)<...<x n (t)} dt = / h i 1 {x 1 (t)<x 2 (t)<...<x n (t)} dr (t) 

«/ «/ 

&i l{x 1 (t)<x 2 (t)<...<x„(t)} dt. 

In addition, the processes /3j(r(-)), i G [n] are martingales with respect to the filtra- 
tion (J>(t))t>o with quadratic covariation processes given by 



l{Xi(t)<x 2 (t)<...<x„(t)} dt, G [n]. 



This identity can be derived by following the steps in the previous display. From the 
last computation we can conclude, in particular, that after extending the underlying 
probability space if necessary, we can find a Brownian motion j3 = (/3i,/3 2 , . . . ,j3 n ) 
with zero drift vector and diffusion matrix € such that 



A(t(-)) = / l{x 1 (t)<x 2 (t)<...<x n {t)}dp i (t), ie[n] 



o 
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Finally, we have 

M^(")) = / 1 {*i(r(*)M J+ i(r(*))} dA i(^(*))= / 1 {x A r(t )) =X ]+l{T{ t))} (dt-dr(t)) 

L{x J (t)=x J+1 (t)} dt - y 1 {x j( t)=x J+1 (t)} d * = ^ IfXfW^+iC*)} dt, 

for j £ [n-1]. Here the second identity is a consequence of f !2.3j) and the fact that 
the boundary local times Af, j' ^ j, do not charge the set {t : Xj(t) = Xj + i(t)} 
(see the main result in [H]); and the fourth identity follows from the fact that the 
instantaneously reflecting Brownian motion Z does not spend time on the boundary 
of the orthant M™" 1 . All in all, we can now conclude that (X, /3) is a weak solution 

to d2HD. 

Step 2. We now turn to the proof of weak uniqueness. To this end, let (X, W) be 
any weak solution to (12. ip . Define 

a{t) = inf js > : jf l{ Xl (a)<x 2 (a)<...<x n (a)} da = tj , t > 0, 

and set X (■) = X(cr(-)). Using Levy's characterization of Brownian motion one 
verifies that 

n-l 

Xi (t) = Xi (0) + bi t + Wi it) + v i,i L 3 (*) » * > o» 

j=l 

where = (Wi, W2, • • • > W n ) is a Brownian motion with zero drift vector and diffu- 
sion matrix £, and {Lj} -^ n _^ are nondecreasing processes whose points of increase 
are contained in the sets 

{t>0: Xj{t) = X j+1 (t)}, jE[n-l], 
respectively. Moreover, the law of X is uniquely determined by the joint law of 

n 

(2.4) (X 2 (•)- Xi (•), X 3 (•)- X 2 (■),..., JM-)-Xn-i(-)) and ^X(-). 

i=l 

However, by the uniqueness result in [T7] we can identify the first of the latter two 
processes as an instantaneously reflected Brownian motion in the orthant so the 

joint law of that process and its boundary local times is uniquely determined. More- 
over, the second process can be constructed by using the first process, its boundary 
local time processes and an additional independent one-dimensional standard Brow- 
nian motion, so the joint law of the processes in ( 12. 4p is uniquely determined. Thus, 
the law of X is uniquely determined as well. Finally, the law of X is also uniquely 
determined as one can verify that X (•) = X (r (•)), where r is defined as in step 1 
above. 



Step 3. Suppose now that Assumption [TJ does not hold. Then a weak solution of (12. ip 
cannot exist. Indeed, if (X, W) was such a weak solution, we could define the time 
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change a(-) as in step 2 above and let X(-) = as before. Then, the arguments 

in step 2 would show that the process of spacings 

(X 2 (•) - X, (■) , X 3 (■) - X 2 (■),..., X n (•) - X n _ x (•) ) 

is a reflecting Brownian motion in the orthant in the sense of [17]. However, 

by the main result in [IT] the latter process does not exist if the reflection matrix Q 
is not completely-iS. This is the desired contradiction. □ 

The following example shows that, even when Assumption [TJ holds, one cannot 
expect a strong solution to (12.11) to exist. 



Example 1. Consider the following specification of parameters: n = 2, b\ = b 2 = 0, 

2 , ^2,1 = \] 



Ci i = C2,2 = 1, Ci 2 = t2,i = 0, V\^i = —5, ^2,1 = 4; that is, the system of SDEs is 



(2-5) dXx(t) = l { x 1 (t)<x 2 (t)} dWi(t) - - l {Xl (t)=x 2 (t)} dt, 

(2.6) dX 2 {t) = l {Xl (t)<x 2 (t)} dW 2 (t) + - l {Xl{ t)=x 2 (t)} dt, 

with Wi, W 2 being independent one-dimensional standard Brownian motions. We 
claim that this system does not admit a strong solution. It is well-known (see The- 
orem 3.2 in [1]) that strong existence and weak uniqueness together imply pathwise 
uniqueness, so it suffices to show that pathwise uniqueness does not hold for the 
system (I2.5p - fl2.6p . To this end, we consider the SDE 

(2.7) dZ{t) = l { z W >o} d/3 (t) + l { z (t )=o} dt, 

where /3 is a Brownian motion with zero drift and diffusion coefficient 2. The main 
result in [5] shows that pathwise uniqueness does not hold for this equation. Therefore 
it suffices to argue that pathwise uniqueness for the system fl2.5|) - fl2.6p would imply 
pathwise uniqueness for the equation f l2.7p . Indeed, let Z, Z' be two solutions of 
(12. 7p on the same probability space and with respect to the same Brownian motion 
p. Extend the probability space so that it supports an independent Brownian motion 
W with zero drift and diffusion coefficient 2, and define S, S' according to 

dS (t) = l {z( t)>o} dW (t) and dS" (t) = V(t)>o } dW(t). 

Finally, set 

S — Z S + Z ! S' — Z , S' + Z' 

Ai = , A 2 = and A\ = , A 9 = . 

2 2 1 2 2 2 

Then both {X\, X 2 ) and (X^X^) are weak solutions of the system fl2.5|) - fl2.6p with 
respect to the Brownian motion ((W — f3) jl, (W + /3) /2). Therefore if pathwise 
uniqueness did hold for the system f l2.5p - fl2.6p . we would be able to conclude that 
Xi = X[ and X 2 = X' 2 pathwise, and, hence, that Z = Z' pathwise; in other words, 
the solution of (12. 7p would be pathwise unique. This is the desired contradiction. 

2.2. Markov property and invariant measures. Having established that the 
weak solution X of the system (12. ip exists and is unique (see Theorem |2]), we can 
now proceed to study some of its properties. First, we remark that weak existence and 
uniqueness imply that the corresponding martingale problem is well-posed (see, e.g., 
Corollary 4.8 and Corollary 4.9 in chapter 5 of p2]). Therefore, by Theorem 6.2.2 
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in [12], the process X is Markovian. In addition, the relation X(-) = X(r(-)), 
where X is an instantaneously reflecting Brownian motion in the wedge W with a 
nondegenerate diffusion matrix, shows that the process X has the Harris property: 

(2.8) Vx, yeW, r>0: F x (\X (t) - y\ < r for some t > 0) > 0. 
Moreover, the corresponding property is true for the process of spacings 

Z (•) = (X 2 (•) - X 1 (•) , X 3 (•) - x 2 (•),..., X n (•) - X n _x (•)) . 

Thus, Z has a unique invariant distribution provided that it is recurrent or, equiva- 
lently, if 

Z (■) = (X 2 (■) - X, (•) , X 3 (■) - x 2 (•),..., X n (•) - X n _! (•) ) 

is recurrent. By Proposition 2.8 in the dissertion [TO], the latter is the case if and 
only if 

(2.9) Q^ 1 (b 2 - bx, b 3 - 62, . . . , K - b n ^) T < 

componentwise. Here, the superscript T stands for the transpose of the vector under 
consideration. We summarize our findings in the next proposition. 

Proposition 3. The processes X , Z are Markovian. Both of them possess the Harris 
property. Moreover, the process Z has a unique invariant distribution if and only if 
the recurrence condition (12. 9p is satisfied. 

For a wide class of coefficients the invariant distribution of the process Z can be 
given explicitly. 

Theorem 4. Suppose that in addition to (12. 9 p the condition 

2A = QD + DQ 

is satisfied, where A is given by (12. 2 j) and D = diag(A) (the diagonal matrix, whose 
diagonal elements coincide with those of A). Then the invariant distribution of the 
process Z is given by 

1 n— 1 

j'=i 

where C = 1+ ?^.' 7 ' +7 i "~ 1 zs i/ie appropriate normalization constant, 

7 = 2 diag (Q) D^ 1 Q^ 1 (b 2 - b u b 3 - b 2 , . . . , b n - b n ^) T , 

I 3 = (7i, • • • > 7i-i, Ti+i> ■ ■ ■ > 7n-i) , i e [n - 1] , 

and dz 3 , j G [n — 1], are the Lebesgue boundary measures on the faces {zj = 0}, 
j £ [n — 1], respectively. 

Proof. Consider the transformed process Z(-) = A~ x l 2 Z(-) which is a sticky Brown- 
ian motion in the cone A -1 / 2 (IR + ) n_1 with drift vector A" 1 ^ 2 (b 2 — b\, . . . ,b n — 6 n _i) T , 
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identity diffusion matrix and reffection matrix A^ 1 ^ 2 Q. Let C be the generator of 

the corresponding diffusion. Then, it is sufficient to check that 

(2.10) 




where v and v\ j G [n — 1], are measures on A 1 / 2 (M+) n 1 and A 1//2 {zj = 0}, j G 

[n — 1], respectively, induced by the measures ^e - ^ 7 '^ dz and "v' ' l{ Zj=0 y dz- 7 , 
j G [n — 1], respectively. From this point on, one just needs to follow the compu- 
tations in the proof of Lemma 2.1 on page 130 in [9| to check that (I2.10p is indeed 
satisfied. □ 



3. Convergence and general setup 

This section is divided into two parts. In the first part (subsection 13. ip we prove 
the convergence theorem (Theorem [1]) as stated in the introduction. Then in the 
second part (subsection 13. 2p we describe a much larger class of particle systems that 
converge to appropriate sticky Brownian motions in W. 



3.1. Proof of the convergence theorem. Given the uniqueness of a weak solution 
to the system of SDEs (I1.3P as proved in Theorem [21 Theorem [1] is a consequence of 
Proposition below. To state and obtain the latter, we study the following decom- 
position. For each i G [n] we can write 
(3.1) 

n— 1 n— 1 n— 1 n—1 

x? 1 a) = x? (o) +a? (t) +y: cff (t) - £ it) +y: Ag M ( t ) - y: &tf ® , 

j=i j=i j=i j=i 



where 



r~iR,M ,.\ ._ -nfl T R,M( + \ ._ nR f i j„ 

^ij W - y y j m W - ^ ^ 1 {xM (s)+ _^ < x« 1 ( s ),xf ( s)+ -^=x« 1 ( s )} as > 



A *j W : / 1 [ X f{s)+^ n <Xf +1 {s),Xf (s) + j-=X^{sj) d ^ 6 Z S ) 



and the processes C^j M , I^f 1 , and A^ M are defined similarly to Cf^ M , lfj M , and 
Ag M , respectively. For m G N, let L> m = D ([0, oo), lR m ). We have the following 
convergence result. 

Proposition 5. Assume that the initial conditions {X M (0) , M > 0} are determin- 
istic and converge to a limit x G W as M — >• oo. T/ien i/je family 

(3.2) { (X M , A M , J L ' M , 7^ M , A L - M , A fl ' M ) , M > 0} 



STICKY BROWNIAN MOTIONS 



11 



is tight in D An2 ~ 2n . Moreover, every limit point (X°°, A°°, I L '°°, I R '°°, A L '°°, A R '°°) 
satisfies the following for each i G [n] and j G [n — 1]: 

n-l 



(3.3) Xf 



1 {xr(-)<...<x«(-)}V^adWi ( s ) + 



3=1 



7 l J "{ x fw =x j+i( s )} ds ' 



(3.4) 
(3.5) 
(3.6) 



tL,00 

i,0 



jR,00 

hi 



1 {x™(s)<X™(s)<-<X™(s)} 



1 {^?°( S )=X- 1 ( S )} ds ' 



2adV^(s] 



A 



L,oo 



A n >: 



R,oo 



0. 



wit/i a suitable n-dimensional standard Brownian motion W = (W\, . . . W n ). 

Proof. Step 1. The tightness of the family in (13 .2p can be verified using the necessary 
and sufficient conditions of Corollary 3.7.4 in [7]. Indeed, note that for i G [n] and 
j E [n — 1], the processes P;(-) - Qi(-), as well as {M 1 ^{R i j{t) - 9 R jt),t> 0) and 
(MV 4 (Ljj-(t) — 9^t),t > 0) all converge to suitable one-dimensional Brownian mo- 
tions in the limit M — > oo. Therefore, the conditions of Corollary 3.7.4 in [7J hold for 
the corresponding families of processes indexed by M > 0. It is now straightforward 
to verify that the same conditions hold for the family (A M I L ' M , I R > M A L ' M \ A R ' M ), 
M > 0, so it is tight on D 4n ~ 3n . In view of the decomposition (13.11) . the first 
statement of the proposition now readily follows. 

Step 2. Now, fix a limit point (X°°, A°°, I L '°°, I R '°°, A L '°°, A R '°°) and to simplify 
notation assume that it is the limit of the whole family (13. 2p as M — > oo. 

We start with a few simple observations about the limit point under consider- 
ation. Note first that, for any fixed M > 0, the jumps of all components of 
(X M , A M , I L ' M , I R ' M , A L,M , A R ' M ) are bounded in absolute value by so all com- 
ponents of the limit point must have continuous paths. Moreover, for every fixed 
t > 0, the family {A M (t) , M > 0} is uniformly integrable due to the estimate 
-t 



E [Af (t) 2 ] = E 



d[Pi-Qi] (s) 



< 2 at, i E\n] 



where [•] denotes the quadratic variation process of a process with paths in D . 
This and the fact that A M is a martingale for any fixed M > show that A°° is a 
martingale with respect to its own filtration. 

Next, we observe that, as limits of non- decreasing processes, and 1^-°° must 
be non- decreasing processes themselves for every i E [n], j E [n — 1]. Finally, for all 
i G [n], j G [n — 1], the quadratic variation processes of the martingales A i ] and 

A R,M ,. r 

A i '• satisfy 



Vt > : 



lim E 

Af->oo 



L,M 
\3 



if) 



and 



lim E 

M->oo 



R,M 



(t) 



0. 
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Therefore, the distributional limits 

A L '°° = lim A L ' M , A R '°° = lim A R ' M 

in D 1 exist and are identically equal to zero. 

Step 3. To show (13.41) . we study the quadratic covariation processes {Xf° 1 X'^) = 
(4°°, Af 3 ), i, i! G [n]. We first claim that (A°°, A°F) = whenever i ^ i'. To this end, 
it suffices to show that for any such pair of indices Af^{-)A'^'{-) is a martingale with 
respect to its own filtration. The latter is the limit in D 1 of the family of martingales 
[Af 1 (-)Af, 1 (■), M > 0} by definition, so it is enough to prove that, for any fixed t > 0, 
the random variables {Af 1 (t)Af?(t), M > 0} are uniformly integrable. The latter is 
a consequence of the following chain of estimates: 



E [Af {tfAf (t) 2 } =E 
= E 
< E 
= E 
= E 



ArUsfdAffis) 2 



Af(s) 2 d[Aff] (s 

f A^isfdiP.-Q,] (s 
Jo 



+ E 


[I 







+ E 


[I 







+ E 


7V(* 


) 2 dAr( S ) 2 




Jo 




\s) 2 d[Af] (s) 




t 







1 V(*) 2 -^ 

Af(s) 2 2ads 

o 



(dp* + dQi,)(s) 
t 



E 



' Af{sf4v7 
o Vm 



(dPi + dQi)(s) 



+ E 


[I 







Af{s) z 2ads 



< 



2a [ E [[Pi - Qi] (s)] + E [[P if - Q it ] (s)] ds < 2a [ 4as ds = 4a 2 t 2 . 
Jo Jo 



We next aim to evaluate {Af), i € [n]. To this end, we note first that, for any fixed 
i E [n] and t > 0, the random variables {Af f (t) 2 ,M > 0} are uniformly integrable 
due to the estimate 



E [4 M (t) 4 ] < E 
+ E 



AA^isfdAf 1 (s 



/ Q A M( s )2^ +4lA M {s)l l_ + _l_ ) diPi + Qi^s) 



Uo 



M 



M M 3 / 2 



The first expectation on the right-hand side is equal to zero, since Af 1 is a martingale 
starting at zero; whereas the second expectation can be bounded above by a constant 
depending only on a and t by arguing as in the preceeding paragraph. Putting this 
together with the fact that the functional 

£2 



(U)1,U 2 , ■ ■ ■ ,U n ) ^ / l{a; 1 (s)<-<w n (s)} ds 
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on D n is lower semicontinuous and using the Portmanteau Theorem, we have that 

„t 2 



E 



G (A°°) (A°° (t 2 ) 2 - 4 00 (t{f - 2, 

G(A°°)(a? {hf-A? (t x f -2a J 
G (A M ) (A? (t 2 ) 2 - Af (h 



limE 

e4-0 



{Xf («)<...<X-»(«)} 

t2 



ds 



l{X-( S )+e<X^_ 1 ( S ), fce[n-l]} ds 



> limsupE 

M->-oo 



*2 



tl 



L {X«( S )+-L<« 1 ( S ),fce[n-l]} 



(1,S 



lim sup E 



'2 



1 {Xt I (s)+-kf<X^ I +1 (s)M[n^l]} d [ P i ~ Q* 



G(A M ) (Af (t 2 f-Af(t{f- 
> limsupE [G (A M ) (Af 1 (t 2 ) 2 - Af 1 (t,) 2 - [Af 1 ] (t 2 ) + [Af 1 ] fa))] = 

M— >oo 

for any nonnegative continuous bounded functional G on D n measurable with respect 
to the a-algebra generated by the coordinate mappings on D n ([0, ii]). Therefore, 
recalling that (X?°) = (Af), we conclude that 

(3.7) V0 < t x < t 2 : (X°°) (t 2 ) - (X?) (h) > 2a [' l {x ~ {s)< ... <x ^ {s)} ds 

holds with probability one. On the other hand, 

E [G (A°°) [AT (t 2 ) 2 - AT it,) 2 - 2a(t 2 - t x ))] 

= lim E [G (A M ) [Af (t 2 ) 2 - Af ( tl ) 2 - 2a(t 2 - tl ))} 

= hm E [G (A M ) (Af (t 2 ) 2 - Af (t,) 2 - - QiKh) + [P t - Qi](h))] < 

for any functional G on D n as above. Hence, 

(3.8) V0 < h < t 2 : (X°°) (t 2 ) - (X?) (ti) < 2a (t 2 - h) 

must hold with probability 1. 

In view of (j3.8|) . we see that in order to improve (13. 7p to an equality, it suffices 
to show that the measure d (X°°) = d (A°°) assigns zero mass to the sets {t > : 
Xj°(t) = X°^ 1 (t)}, j G [n — 1], with probability one. To this end, we first recall that 
for every % G [n] the square integrable martingale Af is the limit in D 1 of the square 
integrable martingales {Af , M > 0}, and the random variables {Af(t) 2 ,M > 0} 
are uniformly integrable for any fixed t > 0. Therefore (A°) is the limit in I* 1 of 
{ [Af] , M > 0}, and so by the Portmanteau Theorem 

E [«A°°) (0 - <A~> W) 2 ] < liminf E [([Af ] (t) - [Af] (t)) 2 

1 



< liminf E 

M->oo 



liminf E 

M->-oo 



1 



(dP t + dQi - dP, - dQ V ) (s] 



j 1 {Xf( S )+-i=<X« 1 ( S ) )te [n-l]} TT d [Pi + Q 4 - P' - Qi/] (s) 
.*/ 



< liminf — -— 



E[Pi(t) + Q;(t) + Pi>(t) + Qi,(t)] = 
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for any fixed i, i' £ [n] and t > with probability one. In view of the path continuity 
of the processes (X°°), i £ [n], this implies 

W> = W> = • • • = (*n) 

with probability one. To conclude the argument, note that the occupation time 
formula for continuous semimartingales (see, e.g., [TS1 Theorem VI. 1.6]) implies that 
the measure 

d (X™! - X°°) = d (X^) + d (Xf ) = 2 d (Xf) = 2 d (Xf°) 

assigns zero mass to the set {t > : X°°(t) = with probability one. 

Hence, equality must hold in (\3.7\i . The representation H3.4[) with a suitable standard 
Brownian motion W = (W\, W 2 , ■ ■ ■ , W n ) now readily follows from the Martingale 
Representation Theorem in the form of Theorem 4.2 in chapter 3 of |12j . 



Step 4- We now turn to the proof of (13.31) . (13. 5p and (13.61) . To this end, for any 
M > 0, i £ [n], j £ [n — 1], we consider the decompositions 

I*f(f) = If \t) - Ig**{t) 

t r t 



iS M (t) = if\t)-i«f>\t) 

It is now easy to check that for each % £ [n], j £ [n — 1], the families of processes 
{lf'\M > 0}, {iJ M ' 2 ,M > 0} and {/g M ' 2 ,M > 0} satisfy the tightness criterion 
of Corollary 3.7.4 in [7]. So, after passing to a subsequence if necessary, we obtain 
the existence of suitable limits /j* 3 ' 1 , -^j°°' 2 and /^j'°°' 2 , respectively; for notational 
convenience we assume that the full families of processes converge jointly to the 
respective limit points. 

Clearly the limits are non- decreasing processes and inherit the properties 

(3.9) V0 < h < t 2 : I-f°*(t 2 ) - ^°°' 2 (ti) < ip\h) - ip'ih), 

(3.10) V0 < h < h : I*f' 2 (t2) - /g°°' 2 (ti) < ip\t 2 ) - ip\h). 

Moreover, arguing as in the second half of the proof of Theorem 4.1 in [20] (see also 
the proof of Proposition 9 in [H]), one obtains 

!{Xf d/°° :1 (t) = 0, 

(3.12) J (l {x ^ 1 ( t )<x^(t)} + ^{x^(t)<x^ +1 (t)}^ dl^°°' 2 (t) = 0, 

(3.13) J (l^ttxxs-^t)} + Mx*>(txx%. 1 W}) d/ 5°°' 2 (*) = 
from the corresponding properties of the prelimit processes. 
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Next, we define the time change 

a(t) = inf js > : J l{Xf>( r )<Xf (r)<...<x«>(r)} dr = t j , t > 

and let X°° (•) = X°° (cr (•)). Using Levy's characterization of Brownian motion, we 
conclude that the components of X°° admit the decomposition 

n—l n— 1 

x°°(-) = xr (o) + m(-) + J2 «* >(•)) + £ (*£■ J 5°°'V(-)) - efe ^rV(-))) 

with W = (Wi, W2, • • • , W n ) being a suitable standard Brownian motion. Moreover, 
arguing as in the proof of Lemma 1 in [IT] (in particular, using the Lyapunov functions 
constructed in the proof of Lemma 4 in [H] as well as fl3.9p - fl3.13p ). we can conclude 
that 

/JrV(-)) = and I*p\a{-)) = 0. 
Furthermore, we can identify the process of spacings 



(i 2 °°(-) - *r (■)> ^(o - ■ ■ ■ , - x 



00 / 

n-l ( 



as a reflected Brownian motion in the orthant (M + ) n_1 with reflection matrix Q 
(defined in Assumption [T]), and the processes /°°' 1 ((t(-)), j G [n— 1], with its boundary 
local times. At this point, to obtain the representations f)3.3p . f)3.5p and (13. 6p . one 
can argue as in step 2 in the proof of Theorem [2j □ 

3.2. General setup. In this last subsection, we introduce a much more general class 
of particle systems which converge to appropriate multidimensional sticky Brownian 
motions in the sense of Theorem [TJ We now allow for non-exponential interarrival 
times between the jumps of the particles and for dependence between the arrival 
times of the jumps for different particles. 

To define this more general class of particle systems, we introduce the following 
parameters: n£N for the number of particles as before; a > 0; Af and Af for % G [n]; 
cff, cff", and cff 1 for i, i' G [n]; and finally 0f • and 0fj for i G [n], j G [n — 1]. We 
fix a value M > of the scaling parameter. The random variables and processes we 
define next all depend on M, but for the sake of readability we mostly do not denote 
this dependence explicitly. 

We let [u L (k), k G N} and {u R (k), k G Kf} be two independent sequences of i.i.d. 
random vectors with values in (0, oo) n (the interarrival times between jumps to the 
left and to the right when there are no collisions), and for i G [n], j G [n — 1], let 
{wfj [k) , k G N} and {wfj (k) , k G N} be two independent families of i.i.d. random 
variables taking values in (0, 00) (the interarrival times between jumps to the left 
and to the right when there is a collision). We assume that 

E[«f(l)]= + > E[«f(l)]= («+^) , 

oov(«f (1)) = c^, cov(nf(l),^(l)) = c^, cov(nf (l),n«(l)) = cff, 

E[<(1)] = (fij,-)" 1 » EW(1)] = (*5) -1 - 
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Next, define the corresponding partial sum processes 

i £ 

£tf(0) = 0, of = 5>f (*) , 0?(o) = o, it? (i) = 



fc=l 



fc=l 



WS(°)=°> (3 = ^(0)=0, W£j (£) — ^2 (k) 

k=l k=l 

for all « 6 [n], j G [n — 1] , and also the corresponding renewal processes 

Sf (t) = max {k > : (jfc) <t}, Sf (t) = max {/c > : (Jfe) < t} , 
7j (t) = max {k>0: (k) < t] , 7g (t) = max {/c > : (fc) < *} 



We now define the particle system for any fixed value of the scaling parameter 
M > according to 



(3.14) 



+ 



j n— 1 

.7=1 



.7 

n-1 



1 >^ 

J=l 



Mt 



for ie [n\. Note that the particle configuration (Xf 1 (t) , X 2 M (t) , . . . , Xjf (t)) is an 
element of the discrete wedge W M for any t > 0. We also remark at this point 



that dSHD generalizes (fPj) . Indeed, if Af = Af = for i G [n], cf^ 



R,R 



whenever i ^ i', = cff 1 = a 1 and cff 1 = for i G [n], and all interarrival 
times above are independent exponential random variables with appropriate means, 
then ( EU) reduces to QT2J). 

For the extension of our convergence theorem to particle systems as in A3. 14)) . we 
need the following moment assumption on the interarrival times between jumps. 

Assumption 2. Assume that there exists 5 > such that 

E \uf (l) 2+5 j \ <oo, 

\2+&] . ttt. f R /i \2+<5 



sup max ( E 

M>0 *e[n] 



«f (1) 2+5 



sup max | E 

M>0*S[n],je[n-l] 



+ E 



< oo. 



Under Assumption [2] we have the following convergence result, which generalizes 
Theorem [TJ 

Theorem 6. Suppose that Assumptions^ and^ hold, and that the initial conditions 
{X M (0) , M > 0} are deterministic and converge to a limit x G W as M — > oo. 
Then the laws of the paths of the particle systems {X M (•) , M > 0} defined in f)3.14p 
converge in D ([0, oo), R") to £/ie /<m> o/ t/ie unique weak solution of the system of 
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SDEs 

(3.15) 



n-1 



dx t (t) = (*)<...<*„(*)} (( A f - A 9 dt + « 3/2 d ^ (*)) + E ^(t)=%i(t)} ^ dt 



3=1 



/or i e [ti], taking values in W and starting from x. Here, W = (W\, W2, ■ ■ ■ , W n ) is 
Brownian motion in MJ 1 with zero drift vector and diffusion matrix given by 



an 



d Vij = Qfj — 6f a is as in Assumption^ 



The existence and uniqueness of a weak solution to the system of SDEs (13. 15)) is 
proven in Theorem HI so Theorem is a consequence of Proposition [7J below, which 
is the appropriate generalization of Proposition in subsection 13.11 



As in subsection 13. 1\ we need to study an appropriate decomposition of the particle 
dynamics. For each i £ [71], we write 



n—1 n—1 n—1 n—1 



X?{t)=X? (0) + Af it) + £ Cff [t] - £ eg" (*) + E Aj A/ (t) - £ A#* (t) 

3=1 3=1 3=1 3=1 

where 

Af ®--7Sl SI W*«w->.*M d W (« •) - s. 1 (« •)) • 

/~<R,M t.\ aR T R,M qR f n j 

W- y i,3 J i,3 -^7 1 {xi'(s) + ^ f <X^ 1 (s),Xf'(s) + ^=X^ 1 (s)} aS 

A ^ M{t) := 7^l 1 { Jf i , 'w+^?< jr ffiW^w+^/fiW} drw (^ s ) ' 
rf 7 -(t) :=iS(t)-efit, 



and the processes Cf^ M , lfj M , A^ M , and T^- are defined similarly to C^ M , //^' M , 

Afj M , and T 4 •, respectively. The following proposition is the appropriate general- 
ization of Proposition [5] to the present framework. 

Proposition 7. Suppose that Assumptions^ and^ hold, and that the initial condi- 
tions {X M (0) ,M > 0} are deterministic and converge to a limit x £ W as M — >■ 00. 
Tnen iae family 

(3.16) { (X M , A M , J L ' M , I*> M , A L ' M , A R - M ) , M > 0} 
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is tight 
satisfies 


in _D 4 ™ 2 2n . Moreover, every limit point [X°° , A°° , I L '°° , 1^°° , 
the following for each i G [n] and j G [n — 1]: 


A L, °°, A R '°°) 


(3.17) 




■)= / 1 K(S) <...<^ W }(W- A 9 d S + a 3 / 2 dW t 

*/ 

,V1 1 

+ ^2 V i,3 1 {x°°(s)=X°l 1 (s)} ds ' 
J=l 1/0 


W) 


(3.18) 


A ( 


\ t tf\R \L\ j 3/2 jtia 

.0 = / 1 {x 1 -(-X...<x-w} (( A * " A J ds + a 7 

«/ 


(*)), 


(3.19) 


tL,00 / 

ij ' 


-) = / 1 { x rw= x i+iW} ds ' 




(3.20) 


tR,oo i 

hi ' 


•) = Ag°°(-) = o, 





a Brownian motion W = (Wi, W 2 , ■ ■ ■ , W n ) as m £/ie statement of Theorem® 

Proof. One can proceed as in the proof of Proposition [5l so we only explain the 
arguments which are different. First, note that Theorem 14.6 in [3] and its proof 
extend to the case of the multidimensional renewal processes 

(<J L \ i^ R \ iT L l (T R \ 

\ u i Jie[n] ' V « Jie[n] ' I «JJi6[n] j6[n-l] ' I M /ie[ra],je[n-l] ' 

yielding the joint convergence of 

{(M-^ftMi) - Sf(Mt)),t > 0)} i6[B] , 

{(M-^T^(VMt),t>0)} 

I V /J ie[n]j'6[n-l] 

and (fM-^Tf (VMt),t > o)) 

IV ,J J ) ie[n] ,j€[n-l] 

to appropriate Brownian motions. The rest of step 1 and step 2 in the proof of 
Proposition [5] carry over to the present setting in a straightforward manner. 

Now, one needs to show that every limit point (X°°, A°°, I L '°°, I R <°°, A L <°°, A R <°°) 
satisfies 

ds. 



(X°°,X™) (•) = (A°°,A™) (•) = a a tV l{ x -( s )<...<x-(s)} 

To this end, one can first proceed as in step 3 in the proof of Proposition [5] to show 
that 

(3.21) d(Xr,X^) = ^d<Xf), i, % G [n], j G [n]. 

Next, one can invoke the Portmanteau Theorem as before to conclude that for all 
i G [n] and < t\ < t 2 \ 

(3.22) a 3 c M f \ {xns)< ... <x ^ s)} ds < (t 2 ) - (X?) (h) < a 3 c M (t 2 - h). 

Jti 
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Moreover, since the measures d(X°°), j G [n — 1], assign zero mass to the sets 
{t > : X°°{t) = X^it)}, j E [n - 1], respectively, f lCTj) and (l3^2|) suffice to 
identify all quadratic covariation processes (Xf°,X°?), G [n\. Similarly, one can 
identify the bounded variation parts of the processes A°°, i G [n], as multiples of the 
quadratic variation processes (X°°), i G [n], respectively. The rest of the proof can 
be carried out by following the arguments in step 4 in the proof of Proposition □ 
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